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Abstract 

We consider a Sinai billiard where the usual hard disk scatterer is replaced by 
a repulsive potential with V{r) ~ Ar~" close to the origin. Using periodic or- 
bit theory and numerical evidence we show that its spectral statistics tends to 
Poisson statistics for large energies when a < 2 and to Wigner-Dyson statis- 
tics when a > 2, while for a = 2 it is independent of energy, but depends 
on A. We apply the approach of Altshuler and Levitov [Phys. Rep. 288, 
487 (1997)] to show that the transition in the spectral statistics is accompa- 
nied by a dynamical localization-delocalization transition. This behaviour is 
reminiscent of a metal- insulator transition in disordered electronic systems. 
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The statistical distribution of quantum mechanical energy eigenvalues is of fundamental 
interest in diverse areas of physics such as condensed matter, atomic, and nuclear physics [IH. 
In the strict quasiclassical limit, where the de Broglie wavelength is much smaller than all 
other lengthscales, the theory of quantum chaos relates the spectral statistics of a quantum 
system to the dynamics of its classical counterpart. While chaotic classical dynamics leads 
to Wigner-Dyson random matrix statistics 0, integrable dynamics generically gives rise 
to Poisson statistics (uncorrelated eigenvalues) [Q. There is overwhelming evidence, both 
experimental and numerical, for these results [Q. The spectral statistics of systems with 
mixed classical dynamics is expected to be described by superposing level sequences with 
Poisson and random matrix statistics, where the respective mean level spacing is determined 
from the size of the corresponding phase space volume with regular or chaotic dynamics . 

However, when the de Broglie wavelength is not the smallest lengthscale in the system, 
the spectral statistics is not solely determined by the classical dynamics. For example, the 
classical motion of an electron in a (three-dimensional) disordered system, such as a metal 
with substantial impurity scattering, can be regarded as completely chaotic. Nevertheless, 
upon variation of the disorder strength at a fixed energy (Fermi energy), the spectral statis- 
tics undergoes a sharp transition from Wigner-Dyson to Poisson statistics 0. Note that 
at the critical point the de Broglie wavelength is of the same order of magnitude as the 
elastic mean free path, a classical lengthscale. The transition in the spectral statistics is ac- 
companied by a transition from extended to localized eigenstates (Anderson metal-insulator 
transition) [0]. In the case of disordered systems one can therefore attribute the deviation 
of the spectral statistics from what one would expect on the basis of the classical dynamics 
to the quantum phenomenon of localization. 

A similar effect on the spectral statistics may be caused by dynamical localization . 
Dynamical localization occurs for example in a circular billiard with a rough boundary, 
where in a certain range of parameters the quantum eigenstates are localized in angular 
momentum space despite that the classically chaotic dynamics leads to diffusive spreading 
of an initial angular momentum distribution [0. As a function of boundary roughness 
there is a crossover between localized and extended eigenstates, which takes place when the 
de Broglie wavelength AdB is of the same order of magnitude as a classical lengthscale set 
by the roughness, {{dR/d9Y) / Rq, where R{9) defines the rough circle in polar coordinates, 
-Ro = (R), and the average is over the angle |]10[. At the same time the level statistics 



changes smoothly from Poisson to Wigner-Dyson statistics |10| , pT| . 

It is the purpose of the present article to show that a dynamical (non-random) system 
may also display a sharp transition in the spectral statistics. This is accompanied by a 
dynamical localization-delocalization transition of the eigenstates and is reminiscent of a 
metal-insulator transition in disordered systems. We thereby extend the list of analogies 
between dynamical and disordered systems [0]. 

We consider a generalization of the well-known Sinai billiard (SB). Sinai proved that the 
free motion of a classical particle being specularly reflected from a disk of radius R inside a 
square with periodic boundary conditions is completely chaotic [0. We modify the SB by 
replacing the disk with the scattering potential 

nr) = |^((^r-0'^<^ , (1) 
^ ^ 1 , r> R ^ ^ 
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where a and A are positive parameters (see also Fig. |l|). Note that in the hmit a — > cxd the 
SB is recovered. In a chaotic bilhard, the replacement of a hard wall by a soft potential 
barrier generically leads to the formation of stable islands |T^, however in the present case 



these are hardly visible in a Poincare surface of section and are not relevant for the effect 
studied here. 

A similar system was considered by Altshuler and Levitov fl^. They focused on the prop- 



erties of the eigenstates and proved the occurence of a dynamical localization-delocaliztion 
transition. In contrast, the emphasis in the present work is on a transition in the level 
statistics and on its semiclassical origin. The transition occurs as a function of a at a = 2, 
irrespective of the value of A. It is a sharp transition in the limit of large energy. Using 
periodic orbit theory we find that it is caused by a competition between regular ('bounc- 
ing ball '-type |T5[) and chaotic orbits. In the second part of this article we show that the 



approach of Altshuler and Levitov may be applied also to the present case. We thereby 
demonstrate that the transition in the spectral statistics is accompanied by a dynamical 
delocalization transition of the eigenfunctions. 

For large energies, the qualitative quantum dynamics of the generalized Sinai billiard 
(GSB) is determined by the interplay of three lengthscales, (i) the de Broglie wavelength 
AdB, (ii) the radius = (E/X + of the classically forbidden area, and (iii) a typical 

length I of an orbit before its direction is randomized by scattering from the potential. 

When a > 2 the potential (|lD effectively acts as a hard wall, since the radial wave 
function then vanishes like exp(— cr near the origin, i.e. 

faster than any power. To estimate the behaviour of the spectral statistics, we approximate 
the potential (|1|) by a hard disk with the energy dependent radius Tc- According to the 



semiclassical theory [0, the structure in the spectrum on the scale SE is determined by 
periodic orbits of length / = hv/6E, where h denotes Planck's constant and v the velocity of 
the particle. For I > I the chaotic orbits dominate and the spectral statistics will be random- 
matrix-like, while for I < I the regular orbits dominate, leading to Poisson-like statistics. 
With e = 5E/A measuring the energy on the scale of the mean level spacing A, we expect 
random-matrix-like behaviour on energy scales e < e and Poisson-like behaviour for e > e, 
where e ~ A/X^b^ and A denotes the area of the billiard. Now / can be estimated by 
/ ~ A/rc, so that we find e ~ Vc/Xds- This ratio scales with energy as Tc/XdB oc 
Consequently, for a < 2 and increasing energy, e tends to zero, while e —* oo for a > 2. 
On this basis we expect a sharp transition between Poisson and random matrix statistics 
at a = 2 in the limit of large energy E. Intuitively, the classically forbidden area becomes 
invisible to quantum mechanics for a < 2, while for a > 2 it becomes more and more 
sizeable. 

To see this explicitly, we apply Berry's periodic orbit theory of bilinear spectral statistics 
||T6| to the SB with energy dependent disk radius Tc. The SB has two types of periodic 
orbits, those that never strike the disk and those that do. The former are marginally stable, 
occur in one-parameter families, and will be referred to as regular orbits, while the latter are 
unstable, isolated, and will be referred to as chaotic orbits. In the following we determine the 
contribution of the regular orbits to the spectral form factor K{t). Denoting the unfolded 
density of states by d{e), the form factor is related to the two-point correlation function 



R{s) = {d{e + ^)d{e-^)) -1 (2) 
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by Fourier transformation, K{t) = ds e^'^**^-R(s). One way to determine the contribu- 
tion of the regular orbits to the density of states is to suitably modify the trace formula 
for the empty billiard. To avoid degeneracies in the spectrum, we use a rectangular instead 
of a quadratic billiard (sidelengths a and b) and quasi-periodic boundary conditions for the 
wavefunction, i.e. ip{x + a,y) = e^'^'"ip{x,y) and analogously for the y-direction |T^. The 
eigenvalues of the rectangle without the disk are Ejk = {2TTj + (j)xY/a'^ + {2nk + (f)y)'^/b'^ (here 
and below /2m = 1). Applying the Poisson summation formula to the density of states 
of the empty rectangle and replacing the resulting Bessel function by its asymptotic form 
leads to a representation as a sum over periodic orbit families. Each family is specified by 
two integers, m and n (positive or negative), denoting the number of traversals across the 
billiard in x- and ?/-direction, respectively. Now including the disk obstructs the path of the 
orbits that violate the condition 2 ImnfciE) < A, where A = ab and Imn = [(ma)'^ + (n6)^]^/^ 
denotes the length of the orbits. The remaining families have to be weighted by the area 
they cover. The oscillatory part of the regular contribution to the unfolded density of states 
is then 

drcgU) = \ll fmn{E) ^ j^"^ U/2_£;l/4 ' (3) 

(■m,n)' j=l ^-^ mn) 

where the first sum is over all primitive orbits and the second over repetitions. In eq. (^ 
we have introduced the actions Smn = ImnVE + nKp^ + n(t)y and the orbit selection function 
fmn{E) = (1 — 2rc{E) Imn/A) 6{A — 2rc{E) Imn), and set the mean level density to its large 
energy limit, so that e = AE/4tt. This formula is a simple generalization of the expression for 
the usual SB given in Ref. [|18|. Substituting Eq. (^ into (0) and retaining only the diagonal 
terms in the double sum over periodic orbits (the diagonal approximation is justified by the 



energy average and known to be exact for regular systems |]T6[) yields the contribution of 
the regular orbits to the form factor. 



(m,,n)' j=l '^'^jl'mnV^ 

This expression may be evaluated approximately by replacing the sums over primitive orbits 
and repetitions by a single sum over all orbits (m, n). In the large E limit, 

Kresir) ~ (1 - cE^rf 6(1 - cE'^t) , (5) 

with c = 2i?A^/". Eq. (|^) shows that the contribution of the regular orbits to the form 
factor tends to zero for a > 2, while for a < 2 the Poisson limit K{t) = 1 is attained. 
The contribution of the chaotic orbits assures that for a > 2 GOE statistics of random 
matrix theory is reached. To see this explicitly would require to go beyond the diagonal 
approximation. 

Next, we verify the above theoretical prediction of a transition in the spectral statistics 
of the GSB numerically. For a < 2, the eigenvalues of the GSB may be calculated by 
diagonalizing the Hamiltonian in the eigenbasis of the empty billiard. However, for a > 2 
the matrix elements diverge and it is natural to work in the eigenbasis of the 1 / r°-potential 
|T9| . In this regime we chose to apply the Korringa-Kohn-Rostoker method, as described in 
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Ref. [|T8| for the case of the ordinary SB. The difference with respect to Ref. |]T8| is that in 
the present case the scattering phase shifts of the potential have to be determined by solving 
the radial Schrodinger equation numerically, while they can be expressed in terms of Bessel 
functions for the usual SB. This makes the present case computationally more demanding. 

Fig. shows the cumulative spacing distribution N{s) = Jq dx P{x), with P{s) denoting 
the nearest-neighbour spacing distribution, for a = 1 and 4 (we used A = 8 for the former, 
A = 0.01 for the latter, and in both cases R = 1 with a billiard of area 47r). Each curve is 
displayed with two reference curves, N{s) = 1 — for Poisson statistics and the Wigner 
surmise N{s) = 1 — exp(— for the GOE. We observe a clear movement of N{s) towards 
the Poisson curve with increasing energy for a = 1 and towards the Wigner surmise for a = 4. 
Considering the relatively weak energy dependence of the spectral statistics, see Eq. (H), the 
limiting distributions are not expected to be reached within the available spectra. We verified 
that for a = 2 the spectral statistics is independent of energy in accordance with Eq. (|^). 

We turn to the properties of the eigenf unctions of the GSB. To keep the discussion gen- 
eral, we consider the GSB in ci > 2 dimensions. Following Ref. ||l4i we map the Schrodinger 



equation for the GSB on a localization problem. Applying Levitov's criterion of a diverging 



number of resonances |^ we then show that a = 2 corresponds to a critical point associated 
with a delocalization transition. 

Using periodic boundary conditions and a billiard of unit volume, the Schrodinger equa- 
tion in momentum representation takes the form 

Ck + J2 ^k-k'Ck' = ^ Ck , (6) 

where k denotes a site in the reciprocal lattice, Ck the Fourier coefficients of the wavefunction, 
Vk those of the potential, and £'k = -l- Vq. Eq. (H) may be interpreted as the Schrodinger 
equation of a particle on a lattice with on-site energies and hopping amplitudes Vk. 



Considering an eigenfunction with energy E, one finds [|1J] that its Fourier coefficients are 
non- vanishing only in an energy shell E — 6E < i?k < E + 6E, where 6E is proportional to A, 
which is assumed to be much smaller than E. Within this quasi- ((i—l)-dimensional shell, the 
on-site energies ii^k are uniformly distributed quasi-random numbers. Since ^ oc 1/ /c*^"" for 
large k, one may interpret Eq. (M) as an Anderson model with long-range power-law hopping 



within the energy shell [Tj]. A delocalization transition occurs when the mean number of 



resonances per site diverges ||20|. A resonance is defined as a pair of sites k,k' that fulfil 
the condition |\4;_k'| > \Ek — E]^'\- If two sites are in resonance, the eigenstates of the 
corresponding two by two eigenvalue problem have amplitudes Ck,Ck' that are comparable 
in magnitude. Informally, they can then be considered as 'linked' and an eigenstate can 
spread along this link. If the mean number of resonances per site is infinite, there can be no 
localization, irrespective of the hopping strength. 

Noting that for l^_k' 6E the probability that the sites k, k' are in resonance is 
simply \V]i-k'\/6E, the mean number of resonances JV with a fixed site k can be estimated 
by summing the probability over the energy shell, 

A/'= E l^k-k'l - -j= fdkk'-'V, . (7) 
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The sum is estimated by an integration over a flat region of radius \J E (we used Afc = 
5E/2\fE for tfie widtli of the energy shell and dropped an energy independent factor). 
Substituting the asymptotic form of Vk for large k into Eq. we find that the mean 
number of resonances per site diverges in the limit of large energy when a > 2, independent 
of the dimension. In the case that was considered by Altshuler and Levitov [|14| the same 
consideration leads to a critical a of one. For a below the critical value, the eigenstates are 
localized with a power-law tail due to direct hopping, independent of A when c? < 3 (because 
the effective lattice, i.e. the energy shell, is two-dimensional for = 3), and also for small A 
when d > 3. 

In summary, we have shown that the GSB displays a sharp transition in the spectral 
statistics, which is caused by a competition between chaotic and regular orbits. It is accom- 
panied by a dynamical delocalization transition of the eigenstates. The investigation helps 
to clarify the relation between disordered systems and quantum chaos. 

It is a pleasure to acknowledge useful discussions with E. B. Bogomolny, O. Bohigas, 
A. Mirlin, F. v. Oppen, C. Schmit, U. Smilansky and H. A. WeidenmiiUer. 
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FIGURES 



V(r) ~ - 




FIG. 1. Potential of the generalized Sinai billiard. 
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FIG. 2. Cumulative nearest-neighbour spacing distribution for the generalized Sinai billiard. 
The energy levels that were used are indicated below the curves (for a = 4 an average over three 
different aspect ratios a/h at constant area ah = Att was taken). With increasing energy a clear 
movement towards Poisson statistics for a = 1 and towards GOE statistics for a = 4 is visible. 
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